LOCAL EXISTENCE AND UNIQUENESS IN THE LARGEST CRITICAL 
SPACE FOR A SURFACE GROWTH MODEL 
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ABSTRACT. We show the existence and uniqueness of solutions (either local 
or global for small data) for an equation arising in different aspects of surface 
growth. Following the work of Koch and Tataru we consider spaces critical 
with respect to scaling and we prove our results in the largest possible critical 
space such that weak solutions are defined. The uniqueness of global weak so- 
lutions remains unfortunately open, unless the initial conditions are sufficiently 
small. 



1. Introduction 

The analysis of mathematical models for the study of surface growth has at- 
tracted a lot of attention in recent years, one can see for example the reviews in 
[1, 11] and numerous recent publications. See for example [22, 23, 21, 7, 18, 17, 
9], which we comment in detail later. 

In this article we consider a model arising in the growth of amorphous sur- 
faces which is described by the following partial differential equation, 

(1.1) 3 t h+ A 2 h + A|VH| 2 = 0. 

on the whole R d or with periodic boundary conditions. The function h(t, •) 
models a height profile at time t > 0, so d = 1 and d = 2 are the physically 
relevant dimensions. In view of this and of Proposition 2.2 we will restrict the 
analysis to the case d ^ 3 throughout this paper (although most of the compu- 
tation holds without restrictions on the dimension). 

Equation (1.1), which is sometimes referred to as a conservative version of the 
Kuramoto-Sivashinsky equation, arises also in several other models for surface 
growth. The two-dimensional version was suggested in [22, 23, 21] as a phe- 
nomenological model for the growth of an amorphous surface (Zr 65 Al 7 5 Cu 2 7,5) 



Date: Cambridge, March 22, 2010. 

2000 Mathematics Subject Classification. 35B33, 35B45, 35B65, 35K55, 35Qxx, 60H15. 

Key words and phrases, surface growth, critical space, uniqueness, regularity. 

This work has been partially supported by the GNAMPA project Studio delle singolarita di al- 
cune equazioni legate a modelli idrodinamici. Part of the work was done at the Newton institute for 
Mathematical Sciences in Cambridge (UK), whose support is gratefully acknowledged, during 
the program "Stochastic partial differential equations". The authors would like to thank Herbert 
Koch for pointing out the approach to the problem. 

1 



2 



D. BLOMKER AND M. ROMITO 



and more recently as a model in surface erosion using ion-beam sputtering 
[7, 18, 17]. The one-dimensional equation appeared as a model for the bound- 
aries of terraces in the epitaxy of Silicon [9]. 

For simplicity of presentation we consider the rescaled version (1.1) with a- 
dimensional length-scales. Furthermore, we have ignored lower order terms 
like the Kuramoto-Sivashinsky term — | VH| 2 or a linear instability given by +AH. 
These terms can easily be incorporated in the result. 

In the physical literature equation (1.1) is usually subject to space-time white 
noise, which we also have neglected for simplicity of presentation. Indeed, us- 
ing the standard method of looking at the difference between h and the stochas- 
tic convolution, the stochastic PDE can be transformed in a random PDE. If the 
stochastic convolution is sufficiently regular, then for each instance of chance 
the path-wise solvability for the stochastic PDE is completely analogous to the 
results presented here and one only needs to consider additional lower order 
terms. This will be done with more details later in Section 5. 

A crucial open problem for equation (1.1), is the fact that the uniqueness of 
global solutions is not known. We remark that numerical experiments do not 
report any problems of blow up, see Hoppe and Nash [14, 13], or the previ- 
ously stated physics literature. Numerical experiments from Blomker, Gugg 
and Raible [4] furthermore indicate a fast convergence of spectral Galerkin meth- 
ods for averaged surface roughness for the stochastic PDE. 

The existence of global weak solutions in dimension d = 1 on bounded do- 
mains has been studied in [4] (see also the references therein), based on spectral 
Galerkin methods. The crucial estimates are energy-type inequalities which al- 
low for uniform bounds on the L 2 -norm. The method has been significantly 
extended by Blomker, Flandoli and Romito [2] in order to verify the existence of 
a solution that defines a Markov process. Winkler and Stein [25] used Rothe's 
method to verify the existence of a global weak solution, this result has been 
recently extended by Winkler [27] to the two-dimensional case, using energy 
type estimates for J e H dx. 

The authors have showed in [5] the uniqueness of local solutions with ini- 
tial values in the critical Hilbert space H 1 / 2 in the one dimensional case. Local 
uniqueness of continuous solutions in W 1,4 for the stochastic PDE in dimension 
d = 1, 2 can be found in [3]. A regularized problem with a cut-off in the nonlin- 
earity in dimension d = 2 has been studied in Hoppe, Linz and Litvinov [12]. 

In this paper we study existence and uniqueness of solutions with initial data 
in a space of BMO-type, which contains all previous spaces where analogous 
results were proved. For periodic boundary conditions it allows for unique local 
solutions with arbitrary initial data in H 1//2 or the space of continuous functions 
C°. 

A weak solution for (1.1) with initial condition h E L 1 1 oc (R d ) is any distribu- 
tion h on R d with locally square integrable gradient Vh e L 2 oc ( [0, oo) x R d ) such 
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that for every smooth and compactly supported test function 4> G C£°([0, oo) x 
R d l, 



(1.2) 



3cJ) 

h(t, x) — — (t, x) dx dt 
Rd dt 

poo 



h(t,x)A 2 ct>(t,x) dxdt + 



o J 



|VH(t,x)| 2 A4)(t,x) dxdt 



R d 



h (x)4)(0, x) dx. 



Note that VH G ^([0, oo)xR d ) implies also H G Lf oc ([0, oo)xR d ) (cf. Lemma2.1) 
and thus all terms in in (1.2) are well-defined, not only in the sense of distribu- 
tions. Moreover, the solution is only defined up to constants. 

Following the remarkable paper by Koch and Tataru [15], this article pro- 
vides a local existence and uniqueness result in the largest critical space, where 
the above stated definition of weak solutions makes sense. As the equation is 
translation invariant (in space) and invariant with respect to the scaling 

(1.3) h(t,x) — >H(A 4 t,Ax), 

we consider the following scaling-aware invariant version of the L 2 oc space for 
the gradient VH, 



(1.4) 



sup 

xeR d ,R>0 







{ R d + 2 t 


. 



Br(x) 



|VH| 2 dy dt 



The paper is organized as follows. In Section 2 we discuss the space defined 
by (1.4) and show an equivalent representation, and its relation to BMO-type 
spaces. Some admissible initial conditions and examples are discussed in Sec- 
tion 3. 

Based on Banach's fixed-point iteration-scheme, Section 4 provides the exis- 
tence and uniqueness results. Section 5 contains some details on the extension 
of such results to the stochastically forced case. We close the paper with Section 
6, where we show smoothness of solutions. 

2. Function spaces 

Recall first the following result, an easy consequence of Poincare's inequality, 
which ensures that all integrals in (1.2) are well defined. 

Lemma 2.1. If u is a distribution on R d such that Vu G L 2 oc ([0,oo) x R d ), then 
u G L 2 OC ([0, 00) x R d ) and thus u G Lf oc ([0, 00), H 1 1 oc (R d )). 

We consider the linear space X° of functions h. with |VH| G L 2 3C ((0, 00) x R d ) 
and thus H G L 1 2 oc ((0, 00) x R d ) such that the quantity ||H|| x o is finite. 

Furthermore, we define the linear space X of functions such that the following 
norm is finite. 



(2.1) 



x 



sup^UVkWHoo}. 



t>0 
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A local in time version of these spaces can be defined for any R > by 

1 



x o := sup 



-d+2 



|Vk(t,y)| 2 dydt 

B r (x) 



x R := sup (t^llVkWHoo). 

for functions k : [0, R 4 ] x R d — > R. Note that we always identify functions that 
differ only by a constant. This is motivated by the fact that the equation is mass- 
conservative, if the total mass J h dx is finite. Furthermore, solutions are only 
defined up to additive constants. 

In order to track the corresponding spaces for initial values, let A = A 2 . Con- 
sider the Green's function G : [0, oo) xR d -}R associated to the operator A, 

where G (t, x) has the Fourier transform (w.r.t. x) G (t, £,) = e - *'^' 4 . By scaling we 
obtain 

G(t,x) =t- d/4 g(xt- 1/4 ), where g(x) = G(l,x). 

The function g is in the Schwartz class since g(£J = e -1 ^' . 

Define the semigroup e tA by the convolution e tA k = G (t, •) * k. Define the 
space 3° of all functions k : R d — > R d such that the bi-caloric extension e~ tA k is 
in X°, endowed with the semi-norm 

||k|| s o:= ||e- tA k|| x o, 

and the space 23 of all functions k : R d — > R d such that 

||k|| s := ||e- tA k|| x 

is finite, endowed with the semi-norm || ■ Define similarly the local versions 
and !B R of these spaces. 
In contrast to the case of Navier-Stokes in dimension three [15], the spaces ¥> 
and CB° (as well as their local counterparts) turn out to be equivalent, as shown 
by the proposition stated below. This makes the analysis simpler than in the 
Navier-Stokes case. This proposition is the only reason, why we restrict to di- 
mension d ^ 3, as we rely for simplicity on the equivalence of B and SO. 

Proposition 2.2. Assume d < 3. Then there are constants C\, c 2 > such that 

(2-2) Ci|| • ||s < || • ||s° < c 2 || ■ ||s 

and 



(2-3) ci|| • ||s R < || • ||s° < c 2 



for every R > 0. 

Moreover, for every d ^ 1, there exists c 3 > such that 

||k||s° ^ c 3||k||BMO(R d )) 

so BMO(R d ) C T>° and in particular L°°(R d ) c B°. 
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For a definition of the space BMO(R d ) of functions of bounded mean oscil- 
lation and its properties, we refer to Stein [24]. Here we only use an equivalent 
norm on BMO(R d ) given by the Carleson measure characterization (see (2.4)). 

Proof. We start by proving (2.2). The inequality on the right holds in any di- 
mension d ^ 1 since it is straightforward to check that there is c > such that 
II ■ || x () ^ c 1 1 • || x- For the inequality on the left, we need to show that for k e 23°, 
x e R d , t > 0, 



t4V(e" tA k)(x)| =t? 



R d 



G(t,x-u)k(y) dy 



^cllkl 



By scaling and translations invariance, it is sufficient to show the statement for 
t = 1 and x = 0. Since 

e -(l-s)A e -sA d 



e- A = 







it follows by the Cauchy Schwartz inequality that 



|V(e- A k)(0)| 



pi 



o 



pi 

( 



pi 



nez d 

neZ d 



< c||k||so y~ 

neZ d 



R d 

d 



g(y(l-s)-5) V (e- sA k)(-y)duds 
g(u(l-s)-^)V(e- sA k)(u)duds 

g(u(l -s) _ 4") 2 du ds 

|V(e- sA k)(y)| 2 dyds 



Bn 



Pi 



Pi 



g(ys ^) 2 duds 



B n 



where B n are the balls of centre 2d 1/f2 n and radius 1 (so that their union covers 
R d ). By a change of variables, 



pi 



I, 



d 

s 2 



g(ys 3) 2 dyds 



pi 



B n 



d 

s 4 



S-V4 Bt 



g(z) 2 dzds. 



First, |IJ ^ C for all n e Z d , as d 3 and g e L 2 (R d ). Note that d ^ 3 is 
necessary, as for e B n we have J s -i/4 Bn g(z) 2 dz | Hgll^ for s | 0. 

For the convergence of the series consider for s e (0, 1) and ^ B n (i.e. 2|n| > 
V&) that 



g(z) 2 dz^ 



R d 



|g(z)|dz-sup{|g(z)| : |z| e S - 1/4 B n } 



^Csup{|g(z)| : |z| > 2|n|d~ 1/2 - 1} 
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which can be bounded by a summable term, since g is in the Schwartz class. 
The inequality (2.3) for the local spaces proceeds similarly. 

Let 4> be in the Schwartz class, with 4) > 0, and set <t> t (x) = t~ d <fy{j). By 
the Carleson measure characterization of BMO (see Theorem 3, Section 4.3 of 
Stein [24]), we have that (up to a constant) 



(2.4) 



Ikll 2 

I MIbmo 



sup 

c£R d ,R>0 









. 



BrM 



-|((Vcf)) t *k)(y)| 2 dudt 



Note that this is an equivalent norm to the standard definition. We could relax 
the conditions on 4), but we are going to use <p = g which satisfies the stronger 
condition g > 0. On the other hand the definition of 23° given above can be 
restated (up to a constant) as 



£0 = SUp 

xeR d ,R>0 



1 



R d+2 



t|((V4>) t *k)(y)| 2 dydt 



with 4> = g, and so ||k|| 2 o ^ C3||k|| BM o(R D )/ as t 2 ^ R 2 . 

3. Examples 



□ 



In view of Theorem 4.3 we wish to discuss for which initial conditions it is 
possible to find R such that the initial condition is small in the 23 R norm. To this 
aim define 

Z = {k : R d -> R : ||1c|| 3r -»■ as R I 0}. 

We later see in Theorem 4.3 that functions in Z correspond to initial conditions, 
where it is possible to solve the equation (1.1) locally for a small time interval. 

The next lemma shows that the whole L°°(R d ), although it is contained in 
23°, is not contained in Z. We will see later that this implies that our method of 
proof fails to provide local uniqueness of solutions for some initial conditions 
in L°°(R d ), although L°°(R d ) c 23 . 

Lemma 3.1. The following statements hold, 

■ there are functions in L°°(R d ) not belonging to Z, 

■ ifk : R d — > R is bounded and uniformly continuous, then k 6 Z, 

■ ifk : R d — > R has bounded gradient on R d , then k e Z. 

Proof. We prove the first statement. Since the Green's function tensorises, it is 
enough to find a counterexample in dimension d = 1. Let k(x) = l[-i,i] [x], then 

t^|(e- tA k) x (x)|= 4[ g'{^Mv)dy = |g(f±i) - g(f^)|, 

hence (choosing x = — 1 and t = R), since g (x) — > as x — > oo, 

't^|(e- tA k) x (x)|) ^ |g(0) - 9(^)1 — »- g(0) > 0. 



sup 
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Assume now that k is bounded and uniformly continuous and fix e > 0. By 
uniform continuity there is 6 > such that |k(x) — k(y)| ^ e for all x, y G R d 
with |x — y\ ^ 6. Since the integral of Vg is zero, 



t 4 |V(e tA k)(x)| 



R d 



t3|z|^6 



Vg(z)(k(x-zt*) -k(x)) dz 
Vg(z)(k(x-zt') -k(x)) dz 

Vg (z) (k(x - zt* ) - k(x) ) dz 



^ 2||k| 



ti|z|<6 

|Vg(z)| dz+ e||Vg|| L i (a d), 



Jt3|z|^6 

hence limsup R ^ ||k|| BR ^ e||Vg|| L i and as e 4- 0, the claim follows. 
Finally let k : R d — > Rbe such that || Vk|| L ~ < oo, then 



t*|V(e tA k)(x)| 



d-l 
t ~ 



g(^)Vk(y) dy 



^ t 4 ||g|| L i||Vk|| L oo, 



andhence ||k|| SR = R 4 ||g|| L i||Vk|| L °o. □ 

Example 3.2. For d = 1 consider k(x) = log |x|. This function has paramount im- 
portance since it is a stationary solution for problem (1.1) (see [5]). It is an inter- 
esting fact that k is neither a weak nor a mild solution (for instance due to The- 
orem 6.1 which ensures smoothness of solutions according to Definition 4.1). 
Here we will show that k e ¥> but k ^ Z. 

Indeed, consider first, 



t* e" 



-tA 



k) x (x) 



R 



g'(z) log |x — t 4 z| dz 



g'(z)(log(t 4 ) +log|t 4 x — z|) dz 



R 



R 



g'(z)log|t 4 x — z| dz 



where we used that the integral over g' is zero. Now substitute x 
obtain 



t 4 x to 



sup{t*|(a x e _tA lc)(x)|) =sup{ g'(z)log|x-z|dz|) = ||9 x e- A k|| c 



sup 



So it is easy to see that ||k||s is finite, but ||k|| 2R is independent of R and does 
not converge to 0. 

Example 3.3. Consider d = 1 and, for <x > 0, k a (x) = |x| a . Then k a ^ ¥>, but 
HKxIIer ~ * 0. So Z may contain certain unbounded functions, which are not in 
23. 
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As in the previous example, 



t*|(9 x e- 



-tA 



X 



g'(z)|x- t3z| a dz 



g'[z)\t-*x-z\*dz 



„ = tf 113 fe~ A k ,, 



, hence ||k a || s — >■ for R — » but 



Thus ||t33 x (e- tA k 

ll^alls = OO. 

Next lemma, together with the main Theorem 4.3, shows that problem (1.1) 
has locally a unique solution for any H d / 2 (R d ) initial conditions. This recovers 
and extends a result proved in dimension d = 1 in [5]. 

Lemma 3.4. The homogeneous space H d/2 (R d ) is contained in Z, where 



H d / 2 (R d ) = |k : R d -> R : ||k|| A 



d 2 



R d 



|£j d |k(£,)| 2 d£,<oo 



and k denotes Fourier transform o/k 

Proof. If t > and x e R d , by using the properties of Fourier transform and 
convolution, 



t*|V(e tA k)(x)| 



1 

tf J 



t* 



Rd 



Rd 



(Vg)((x-y)-H)k(y)dy 



Rd 



(Vg^x-y^k^e 1 ^ dy dE, 



R d J 



R'l 



Vg(y)e 



k(£,)e it x dy d£ 



R d 



|£,| |Tc(£,)|e" 



l^l 4 t 



df,. 



Given a > 0, split the integral in the last line of formula above in two pieces 
and (R), corresponding to the domains of integration {|£J ^ a} and {|£J > a} 
respectively. We estimate using the Cauchy-Schwarz inequality, 



0<t* 



IfcMO^fel) d^^t5||k|| Hd/2 



|^| 2 - d d£,)^ ^cat4||k|| H 



d/2, 



while by a change of variables and Cauchy-Schwarz' inequality again, 



< t* v 

|f,l>a 

In conclusion 



\t\ d \m)\ 2 dt 



2-d Q -2t|£,| 4 



Br ^ caR 



lHd/2 



df, 



l^l d |k(£,)| 2 d^ 



+ C 



|£J d |k(£,)| 2 d£ 



so we see that limsup R ^ ||k|| 2R is bounded by a quantity which converges to 
as a | oo. □ 



largest critical space for surface growth 
4. The fixed point argument 



Define the map 



r(h,k)(t) = 



,-(t-s)A 



Vh(s) • Vk(s)) ds 



and set 
(4.1) 



J?(H)(t) =e- t/K ho-y{h,h){t). 



We will use the following concept of a mild solution, which is given as a solution 
of the variation of constants formula in (4.2). 

Definition 4.1. We say that h G X solves (1.1) with initial condition ho G S, if 
for all t > 

(4.2) h(t) = e~ tA h - Y[h, h) ft) . 

We call h G X R a local solution, if (4.2) holds only for t G [0, R 4 ]. 

The following Lemma is crucial for the proof of uniqueness and existence. It 
verifies that the nonlinear part is locally Lipschitz. 

Lemma 4.2. The map Y is bi-linear continuous from X x X to X and from X R x X R 
to X R ,for all R > 0. 

Proof The bilinearity is obvious. For the boundedness let x G R d and t > 0, 
then 



ivr(h,k)(t,x) 



J 



VAG(t-s,x-u)Vh(s,-y)Vk(s,-y) dy ds 



ft - S 4 



R' 1 



(VAg) 



x-y 



^ IINIx||k|| x 



o ft - sI^a/s J 



t-s) 1 ^ 
(VAg) 



Vh(s,u)Vkfs,-y)duds 
x-y 



:t-s)V4 



dy ds 



< ||h||x||lc||x||g||w3.i(Rd) 



1 



o (t-s)V*Vs 



ds 



^ t 1/4 Bfi, T)IIH||x||ic||Dc||g||w3.i(Rd ) 



■2' 4' 



= C4L 



-1/4 



where B is the Beta function. The corresponding inequality for the local space 
X R proceeds similarly. □ 

Using the previous Lemma, we can now state and prove our main result. The 
first part states global existence of unique solutions, while the second part is 
about local existence of solutions. 
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Theorem 4.3. There is 8 > such that if ||H || s o ^ 5, then there exists a unique 
(global) solution in X of (1.1) with initial condition Hq. 

Moreover, if ||H || s o ^ 8, then there is a unique local solution in X R of (1.1) on 
[0, R 4 ] with initial condition h . 

Finally, /J ho is periodic and small in S^/or some R > (or it is small in then 
the solution is also periodic. 

In particular, ||ho|| 3 o ^ 5 is true for a suitable value of R for all h G Z. 

Proof. We prove the first statement by a fixed point iteration argument. Let c 4 
be the constant defined in the proof of Lemma 4.2 and choose 6 > 0, K > such 
that 

4c 4 S 1 / / 4c 4 6\ „ 1 

1 — >0, 1-A/l — )^K< . 

ci 2c 4 V V Ci / 2c 4 

Define 

(4.3) H = 0, H n+1 = J?(H n ) = e- tA h - Y(H n , HJ, 

then || Hi || x ^ and it is easy to check by induction (and by the choice of 6 and 
K) that ||H n || x 4 K for all n. Then 

||H n+ i — H n ||x = ||^"(H n , H n ) — y(H n _i, H n _i) ||x ^ 2c 4 K||H n — H u _i||x 

and so (H n ) neN is convergent in X to a fixed point of & . 

The same proof works for local spaces, since both constants c.\ and c 4 do not 
depend on R. Finally, if ho is periodic, the statement follows by translation 
invariance and uniqueness. □ 

Remark 4.4 (Forward self-similar solutions). The theorem above allows to show 
the existence of self-similar solutions, namely solutions invariant for the scal- 
ing (1.3). Indeed, assume to have ho G B° (or in a local space) such that h (Ax) = 
ho (x) for all A > 0, then it is easy to verify that H 4 is invariant for the scaling (1.3) 
and that h) is also invariant if so is h. In conclusion the whole sequence 
(H n ) neN defined in (4.3) is invariant, as well as its limit. 

Given a (forward) self-similar solution h, one can write h(t,x) = ^(x/t 1 / 4 ), 
where 4>(x) = h(l, x) solves the equation 

A 2 -4> + A|V\b| 2 - ^x • = 0. 

The simplest case corresponds to d = 1, where the only admissible initial con- 
ditions are all functions ho constants on (— oo, 0) and on (0, oo) (possibly with 
different values on the two half-lines), with ||h ||s R = |ho(l) — ho (—1)1 ||g||i_°°- 

Backward self-similar solutions might provide examples of solutions with 
blow-up. Due to the scaling of the problem, the quantity blowing up is related 
to the derivative of the solution. We do not know if backward self-similar so- 
lutions exist (notice that backward self-similar solutions do not exist for the 
Navier-Stokes equations, see [19]). 
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5. THE STOCHASTIC PROBLEM 

In this section we give a short outline of the proof of local existence for the 
stochastic PDE, without many details on probability theory For details we refer 
to [8, 6, 16]. Consider 

(5.1) 3 t h + A 2 h + A|Vh| 2 = 3 t W, 



where 3 t W is the generalized derivative of a Hilbert-space value Wiener pro- 
cess. Define the corresponding Ornstein-Uhlenbeck process for t > as the 
following Ito-integral 



(5.2) 



Zft) 



,-ft-s)A 



dW 



Note that Z solves 9 t Z + A 2 Z = 9 t W with Z(0) = 0. The mild solution of (5.1) 
is analogous to Definition 4.1 given by a solution of 



h(t) = e~ tA h - y (h, h) (t) + Z(t) 



Now the main problem in the stochastic setting is to determine the regularity 
of Z. Once we know this, we can solve the equation using Banach's fixed point 
argument, as in Theorem 4.3. Moreover v = h — Z solves the following random 
PDE 

(5.3) 3 t v + A 2 v + A| Vv| 2 = -A|VZ| 2 - 2A|Vv • VZ| 2 , v(0)=ho, 

which only contains lower order terms that do not change the proofs, once Z is 
sufficiently regular. 

In the case of bounded intervals (i. e. d = 1) with periodic boundary condi- 
tions and space-time white noise the stochastic convolution Z and its deriva- 
tive 3 X Z are continuous in both space and time, which can be verified using the 
methods in [8]. See for example [20]. This implies that almost surely ||Z||x R — > 
for R — > 0, and we can solve the stochastic PDE (5.2) (or the random PDE (5.3)) 
uniquely in Xr, for some small (random) R > if the initial condition h is such 
that the PDE (1.1) has a unique local solution. 

An interesting question appears in the case of periodic boundary conditions 
and d = 2, as for space-time white noise the convolution Z just fails to be differ- 
entiable in space. Nevertheless, Z will be differentiable, if we consider slightly 
more regular noise. 

For stochastic PDEs on unbounded domains one can use the formulation of 
Walsh [26], although one has to consider that for space-time white noise the 
stochastic convolution Z(t, x) is unbounded for |x| — > oo. 
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6. Smoothness of solutions 

Following the same methods of [10], we show that solutions in X (or X R ) are 
smooth. Define for m ^ 1, 

||k|| x ,m :=sup{t^ V HD^Iloo) 
t>0 L i I — , i J 

and denote by 1 1 • 1 1 x R , m the corresponding local version, where for <x = ( cti , . . . , <Xd) 

we used D a = 3*^ . . . 3** and |a| = <Xi + h a d . 

Let X m be the space 

X m = {k : R d -> R : ||k|| X m := max ||k|| x j < oo}, 

and denote by Xj^ 1 the corresponding local version. For simplicity of notations 
we understand that || • ||x,o = II • II x and for R = oo that X™ = X m . The main the- 
orem of this section is the following result on smoothness in space. Smoothness 
in time then follows from the PDE by a standard bootstrapping argument. 

Theorem 6.1. Let h be a solution of (1.1) in X R , with < R ^ oo. Then h[t) e 
C^{R d ) for all t e (0,R). 

Proof If the initial condition is small enough in B^, the statement follows from 
Proposition 6.2 below. In the general case we notice that if h e Xr, then VH(t) 
is bounded for all t e (0, R), therefore h(t) e Z, by Lemma 3.1. The conclusion 
then follows again from Proposition 6.2. □ 

In order to complete the proof of the above theorem, we need the following 
proposition, which gives also a better estimate of the solution near t = if the 
initial condition is small enough. 

Proposition 6.2. There exists 8 > such that if||rLo|| s o < 5, then the solution to (1.1) 
granted by Theorem 4.3 is in X m for all m ^ 1. 

IfR > and ||H || s o < 5, then the solution to (1.1) granted by Theorem 4.3 is in 
for all m ^ 1. 

We start by giving a slight generalization of (2.2) and (2.3). 
Lemma 6.3. Let < R ^ oo and k e 25 R , then for every m ^ 0, 

(6.1) sup{t^ Y_ ll Da ^ tAk )L} ^cm d (m+l)^||Vg||^ (Rd) ||k|| s o. 

t: ^ R |a|=m+l 

Proof. Since for |<x| = m + 1, 

d 

D a (e~ tA k) = n^ 9 "i e-^ tA )k, 

i=l 
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it is sufficient to show that the operator 3 X . e~ tA maps L°°(R d ) into itself with 
operator norm ||3 X . e~ tA ||l°°^l«> ^ t~ 1/4 || Vg|| L i( R dj. This is immediate since by 
a change of variables, 



tt|3 Xi ( e - tA k)(x)| = tt 



;a xi g)(^f)k(y)dy ^ ||Vg|| L i (a d; 



Finally, #({cc : \oc\ = m + 1}) = ^ cm d . 

Lemma 6.4. There is c 5 > such that for m^l, 0<R^oo and h,ke X£\ 



□ 



||r(H,k)||x R ,m ^ c 5 m d (m+ l) J ^||Vg||P (Rd) ||g|| W 3,i(Rd ) ||h||x R 



S R 



(6.2) 



C5||h||x R ||lc||x R , m +c 5 ||H|| x ||k|| XR 



m-1 



X R J II MIX R ,m-j • 



Proof. Fix m^l, < R ^ oo, t^R and h, k e X£\ Consider a value e e (0, 1) 
which will be specified later, and let |cc| = m + 1. Since |a| ^ 1, there is i ^ d 
such that cti > 1. So assume without loss of generality that ai ^ 1 and let 

a' = a- (1,0, ... ,0). 



D a r(h, k)(t) 



a a /^-(t-slA 



D a A(e 



pt(l-e) 



(Vh(s)Vk(s))) ds 



D a A( e - (t - s]A (Vh(s)Vk(s))) ds 



t(l-e) 



D a A( e - (t - s)A (VH(s)Vk(s))) ds 



For the term (T) we use the factorization introduced in the previous lemma and 
we proceed as in the proof of Lemma 4.2, 



m+3 



ft(l-e) 



®K (m+l)— ||Vg||^||3 Xl Ag|| Ll ||h||x R ||k||x R 



m + 3 



^2v / T^I(^ i ) 4 t-^||Vg||p||3 Xl Ag|| L i||h| 
For the second term we use Leibniz formula, 



- 1 / m+3 , 

s 2 (t — s) 4 ds 



k| 



X R MX 



ft 



t(l-e) 



3 Xl Ae- (t - s)A (D p VH)(D a '- p Vk) ds 
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and, as in the proof of Lemma 4.2, 



4e^ m+ i ^— / <x'\ „, , 

< (1 _ e) ^ t-r- HB^Aglk, Z ( p J INk,,,, II* 



r * ! ds 

m+2 , , 3 US 

t(l-e) S 4 (t — SJ4 



X 



R,m-|PI 



If we set e = ( m+ 1 d) 4d the term 4e 1/4 (l—e) (m+2 ^ /4 (ra+d) d is uniformly bounded 
in m (we recall that the number of multi-indices cc such that |a| = m + 1 is 
bounded by (ra + d) d ) and so by summing up over a the estimates for and 
(2) together show the lemma. □ 

As in the proof of Theorem 4.3, define H = and 

H n+1 (t) = e- tA ho - nH n , H n )(t). 

Lemma 6.5. There is 5' > such that ifO < R < oo and ||h || 2 o < 5, then for every 
m ^ there is K m > such that 

||H n ||x R ,m ^ Km- 

Let us remark that with an explicit estimate of the constants K m , and in par- 
ticular their growth in terms of ra, one could show that solutions are analytic in 
space. For simplicity of presentation, we will not focus on this. 

Proof. If ||ho||so is small enough, the proof of Theorem 4.3 shows that there 
is K such that ||H n ||x R ^ K . By possibly taking ||h || 2 o smaller, we can as- 
sume that A = 2c 5 K < 1, where c 5 is given in Lemma 6.4. We prove the 
statement by induction: the case m = has been already proved. Set a m = 
c 5 m d (m + l) (m+3)/2 ||Vg||p:||g||w 3 . 1 (this is the coefficient appearing in the first 
line of formula (6.2)) and b m = cm d (m + i)( m+1 )/4|| Vg||™ (this appears in for- 
mula (6.1)), then by Lemmas 6.3 and 6.4, 

l|H n + 1 || X R ,1Tl ^ ||Hi||x R ,m+ ||^(H n , H n )|| 

X R ,m 

^ b m K + a m KQ + 2c5Ko||H n ||x R ,m 

m-1 



+ c 5 m d Y_ (7)H Hn ll 3CrJ II II X R .m— j 
5=1 

m-1 

^ (b m K + a m K 2 + c 5 m d (^KjK^j) + A||H 



ti Ii.mi 
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so that by recurrence and again Lemma 6.3, 

||H n+ i||x R ,m ^ 

m-1 

*C (b m K + a m K 2 + c 5 m d (?) K j K m-j) (1 + • ■ ■ + A"" 1 ) + A n ||Hi||x R , m 

5=1 

1 m— 1 

^ T^A ( bmK ° + amK o + Csmd T- (T) K i K ^-j) + b ^ K o> 

j=l 

and the last line in the formula above provides K m . □ 

Proof of Proposition 6.2. Theorem 4.3 ensures that if ||ho||s Q is small enough, then 
there is A = 2c 4 K < 1 (where the number K is given by previous lemma) such 
that ||H n+ i — H n ||x R ^ cA™. We prove by induction that there are numbers 
C m > and \i e (0, 1) such that 

||H n+ i — H n ||x R , m ^ C m |x n , m ^ 0, 

if ||rio||s° is small enough. Let A = 2c 5 K (where c 5 has been introduced in 
Lemma 6.4), assume A < 1 and let A < \i< 1. We have already verified that the 
inductive claim is true for m = 0. Assume the claim is true for 0, . . . , m — 1, then 
by Lemma 6.4 and the inductive assumption, 

ll^n+l — H n ||x R! m 

^ ||^(H n , H u — H n _i)||x R ,m + || ^(H n — Htv-i, H n 

— 1 ) || X R ,m 

< [a m (||H n || XR + ||H n -i||x R ) + C 5 (||H n ||x R ,m + ||H n _i||x R ,m)] ||H n — H n _i||x R 
+ c 5(||Hn||x R + ||H n _i||x R )||H n — H n _ 1 ||x R ,m 
m-1 

+ c 5 m d ^ ^(llH^llx^ + IIH^Hx^OllH^-Hn-illx^H 

j=i 

^ A||H n — H n _i||x R ,m + KmM- n , 

where we have set a m = c 5 m d (m + l) (m+3:i/2 ||Vg||p:||g||w 3 . 1 (the coefficient in 
the first line of (6.2)), K m = 2C (a m K + c 5 K m ) + 2c 5 m d (?)KjC mH/ and 

the constants Kj are given by the previous lemma. By recurrence (notice that 
\x > A), it is easy to see that for every n, 

||H n+1 - H n ||x R , m < A n+1 ||H! - H ||x R ,m + Km(A n_1 + A n ~V + • ■ ■ + ^) 

< (AK m + ^K m )^t n , 

which concludes the induction. In conclusion, the sequence (H n ) neN converges 
in all spaces X£\ □ 
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